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ABSTRACT 

The  sensitivities  of  the  longwave  radiation  emitted  by  the  Earth- 
atmosphere  system  to  variations  in  surface  temperature  and  cloud  amount  are 
evaluated  from  observations.   The  technique  used  is  linear  regression  applied 
to  zonal  average  data  sets  of  longwave  radiation,  cloud  amounts,  and  surface 
temperature.   The  dependence  of  these  sensitivity  coefficients  on  different 
radiation  and  cloud  data  sets  and  on  the  time  averaging  period  is  examined. 
The  results  show  a  strong  dependence  of  the  sensitivity  coefficient  to  the 
specific  radiation  or  cloud  amount  data  set  used  in  the  analysis  and  also  on 
whether  one  uses  annual  or  monthly  means.   To  a  much  smaller  extent  this  is 
also  true  of  the  temperature  sensitivity  coefficient. 

Based  on  simulations,  a  simple  linear  relationship  is  derived  between  the 
planetary  albedo  and  the  surface  albedo  for  the  case  of  clear  skies.   This 
relationship  is  checked  by  applying  it  to  zonally  averaged  surface  albedos  in 
order  to  predict  the  clear-sky  planetary  albedos,  which  are  then  compared  to 
satellite  determinations  of  zonally  averaged  minimum  albedos.   The  results 
indicate  that  this  simple  relationship  can  be  used  to  estimate  clear-sky 
planetary  albedo  to  ±  0.05.   Likewise,  application  of  this  kind  of 
relationship  to  satellite  minimum  albedo  observations  of  the  Sahara  desert 
region  yields  a  surface  albedo  of  0.42,  considerably  higher  than  previous 
estimates. 


T.  S.  Chen,  National  Earth  Satellite  Service,  NOAA,  Washington,  DC  20233; 
George  Ohring,  Department  of  Geophysics  and  Planetary  Sciences,  Tel  Aviv 
University,  Tel  Aviv,  Israel  and  National  Earth  Satellite  Service,  NOAA, 
Washington,  DC  20233;  and  Haim  Ganot,  Department  of  Geophysics  and  Planetary 
Sciences,  Tel  Aviv  University,  Tel  Aviv,  Israel. 


1.    Introduction 

Largely  as  a  result  of  the  cloud-radiation  feedback  problem  in  climate 
studies  there  have  been  a  number  of  recent  attempts  to  use  satellite  observations 
to  evaluate  the  effects  of  clouds  on  the  radiation  budget  of  the  Earth-atmosphere 
system.   Studies  with  some  of  the  broad-band  observations  (e.g.  Cess,  1976;  Cess 
et  al. ,  1982)  indicate  that  for  the  earth  as  a  whole  the  net  effect  of  clouds  is 
negligible,  their  albedo  effect  being  approximately  balanced  by  their  greenhouse 
effect.   Studies  with  radiation  budget  estimates  from  narrow-band  observations 
(Ohring  and  Clapp,  1980;  Hartmann  and  Short,  1980;  Ohring  et  al. ,  1981)  suggest 
that  the  albedo  effect  of  clouds  dominates  their  greenhouse  effect  so  that  the  net 
effect  of  clouds  on  the  radiation  budget  is  not  at  all  negligible.   Both  Hartmann 
and  Short  (1980)  and  Ohring  et  al. ,  (1981)  point  out  the  importance  of  the  large 
geographical  variations  in  the  cloud  effects. 

While  the  issue  of  cloud  impact  on  the  net  radiation  appears  to  be  unsettled, 
there  is  little  doubt  that  the  determination  of  cloud  sensitivity  parameters 
(Schneider,  1972)  depends,  in  large  part,  on  accurate  estimates  of  the  emitted 
longwave  radiation  and  cloud  amounts.   The  sensitivities  of  longwave  radiation 
to  changes  in  cloud  amount  and  to  changes  in  surface  temperature  constitute  two 
of  the  most  important  parameters  used  in  energy  balance  climate  models.   Since 
additional  satellite  data  have  become  available  in  recent  years,  a  further  exami- 
nation of  the  sensitivities  of  emitted  longwave  radiation  to  both  surface  tempera- 
ture and  cloud  seems  readily  justified.   This  paper  addresses  primarily  the  obser- 
vational determination  of  these  sensitivity  parameters,  and  their  dependence  on 
satellite  radiation  data  set.   Some  earlier  work  along  these  lines  has  been  reported 
by  Warren  and  Schneider  (1979) ,  Ohring  and  Clapp  (1980) ,  and  Simmonds  and 
Chidzey  (1982). 
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Another  objective  of  this  study  is  to  develop  and  test  a  simple 
relationship  between  planetary  albedo  and  surface  albedo  for  the  case  of  clear 
skies.   Such  a  relationship  has  a  number  of  important  applications  including 
climate  modeling,  cloudiness  determinations  from  satellite  observations,  and 
surface  energy  budget  estimates  from  satellite  observations. 

2.   Methods  and  data 

Budyko  (1969)  suggested  a  temperature  and  cloud  dependent  longwave 
parameterization  scheme  for  use  in  climate  modeling  studies: 

F  =  a+bTg  -  (aj  +bTTs)  Ac  (1) 

where  a,  b,  a,,  and  b,,  are  constants,  F  is  the  emitted  longwave  flux,  T  is 
the  surface  temperature,  and  A  is  the  cloud  amount.   The  relationship  is 
intended  to  be  applied  to  large  scale  average  data,  e.g.,  to  monthly  means. 
Cess  (1976)  modified  Eq .  (1)  to 

F  =  a+bTs  +  cAc  (2) 

The  sensitivities  of  emitted  longwave  radiation  to  surface  temperature  and 
cloud  amount  are  equal  to  b  and  c  respectively,  and  these  may  be  determined 
from  a  regression  equation  derived  from  data  sets  of  F,  T  ,  and  A  . 

Three  sets  of  longwave  radiation  data  were  employed  in  the  present 
study.   They  are:   2  years  of  NIMBUS  6  data  (Campbell  and  Vonder  Haar,  1980), 
the  pre-1972  satellite  data  (Ellis  and  Vonder  Haar,  1976),  and  the  NOAA  SR 
data  from  June  1974  to  February  1978  (Gruber  and  Winston,  1978).   The  data  are 
zonally  averaged  at  an  interval  of  10°  latitude.  In  addition,  both  NOAA  SR  and 
NIMBUS  6  data  are  also  available  in  zonal  averages  for  land  and  ocean 
separately. 
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The  zonal  monthly  mean  temperatures  are  taken  from  Warren  and  Schneider 
(1979);  the  zonal  monthly  land  and  ocean  mean  temperatures  are  from  Robock 
(1980).   Two  different  sets  of  cloud  data  are  used:   1)  Berlyand  et  al. 
(1980),  which  contains  annual  and  monthly  zonal  means  and  also  includes  values 
for  land  and  ocean  separately,  and  2)  London  (1957),  which  contains  annual, 
zonal  means.   In  deriving  regression  equation  (2)  from  these  data  sets  it  is 
implicitly  assumed  that  each  data  set  represents  the  long  term  mean  for  the 
variable  that  it  represents. 

The  relationship  between  planetary  albedo  and  surface  albedo  for  clear 
sky  conditions  is  developed  from  the  solar  radiation  absorption  model 
described  by  Lacis  and  Hansen  (1974).   In  their  approach,  absorption  of  solar 
radiation  is  parameterized  as  a  function  of  the  water  vapor  distribution,  the 
zenith  angle  of  the  sun,  the  albedo  of  the  earth's  surface  and  the  ozone 
distribution.   With  this  treatment  of  solar  radiation,  the  planetary  albedo 
was  generated  for  a  variety  of  surface  albedos,  solar  zenith  angles  and  water 
vapor  distributions.   Multiple  regression  analysis  revealed  that  the  planetary 
albedo  is  mainly  a  linear  function  of  surface  albedo.   For  example,  based  on 
1080  data  points,  surface  albedo  alone  was  able  to  explain  98%  of  the  variance 
of  planetary  albedo,  with  a  standard  error  of  planetary  albedo  of  0.028.   Thus 
a  one-predictor  regression  scheme  using  surface  albedo  alone  is  justified.   To 
assess  the  reliability  of  the  calculated  relationships,  the  derived  regression 
equation  was  tested  using  observed  data.   They  consist  of  zonally  averaged 
minimum  planetary  albedo  derived  from  NOAA  SR  (assumed  to  be  equal  to  the 
clear  sky  albedo),  zonally  averaged  surface  albedo  by  Robock  (1980),  and  the 
so-called  clear  sky,  local-noon  surface  albedo  by  Hummel  and  Reck  (1979). 
Surface  albedos  from  both  Robock  and  Hummel  and  Reck  were  read  off  from  their 
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figures  for  the  months  of  January,  April,  July  and  October.   In  the  case  of 
Hummel  and  Reek's  data,  their  seasonal  values  were  assumed  to  represent  the 
mid-seasonal  months.   Further,  in  order  to  avoid  certain  uncertainties  in  the 
data  at  higher  latitudes,  all  data  used  are  restricted  to  60°N  to  60°S . 
3.  Discussion  of  results 

3.1  Estimates  of  sensitivity  parameters 

Table  1  lists  the  coefficients  a,  b,  and  c  and  their  associated  statistical 
variables  for  the  annually  averaged  case.   When  temperature  is  considered  the 
only  independent  variable,  the  coefficient  b=(dF/dT  )  ranges  from  1.60  to  1.73  for 
the  global  average,  1.57  to  1.99  for  the  Northern  Hemisphere  and  1.60  to  1.76  for 
the  Southern  Hemisphere  depending  on  the  data  sets.  In  general,  b  shows  less 
spatial  variations  using  NOAA's  SR  and  NIMBUS  6  data.   When  the  fractional  cloud 
cover  is  included,  the  coefficient  b,  which  is  now  a  partial  derivative  (8F/9TS), 
is  only  slightly  less  than  the  value  of  the  total  derivative  dF/dTs  for  the  globe 
as  a  whole  and  for  the  Southern  Hemisphere.   But  for  the  Northern  Hemisphere 
8F/8TC  is  much  less  than  dF/dT  •  this  small  value  of  9F/9T   in  the  Northern  Hemis- 
phere  is  associated  with  a  large  (absolute  magnitude)  value  of  c=(3F/9A  ).   Overall, 
the  values  of  c  vary  substantially  with  longwave  flux  data  set,  suggesting  that 
very  accurate  flux  measurements  are  required.   A  similar  result  was  obtained  by 
Simmonds  and  Chidzey  (1982),  who  examined  the  Ellis  and  Vonder  Haar  (1976)  and 
NOAA-SR  radiation  data  sets  with  the  use  of  the  London  (1957)  cloud  data. 

Listed  in  Table  2  are  results  based  upon  the  same  sets  of  data  except  that 
they  now  pertain  to  monthly  rather  then  annual  averages.   When  they  are  compared 
with  Table  1,  the  values  of  b=(dF/dTs)  appear  to  be  very  consistent,  but  the  values 
of  c  are  considerably  smaller,  particularly  in  the  Northern  Hemisphere  where  a  reduc- 
tion by  nearly  a  factor  of  3  occurs.   These  sensitivities  should  be  more  reliable 


than  those  computed  from  mean  annual  data  since  they  include  the  effect  of 
seasonal  as  well  as  latitudinal  variations.  Similar  qualitative  differences 
between  coefficients  based  on  seasonal  and  annual  data  were  found  by  Simmonds  and 
Chidzey  (1982),  although  they  used  only  the  four  mid-seasonal  months  to  portray 
the  annual  cycle. 

To  further  understand  the  possible  cause  for  the  large  reductions  in  c, 
as  one  goes  from  mean  annual  to  monthly  data,  especially  in  the  Northern  Hemisphere, 
computations  were  made  separately  for  land  and  ocean.   As  seen  from  Table  3,  the 
values  of b=(dF/dT  )  are  again  very  consistent  between  monthly  and  annual  averages 
for  both  land  and  ocean.   But  when  cloud  fraction  is  introduced,  a  large  decrease 
in  c  and  a  relatively  smaller  incerease  in  b=(3F/3Ts)  take  place  simultaneously 
as  one  goes  from  annual  to  monthly  mean  data  in  the  ocean  case.   This  is  in  line 
with  what  we  found  earlier  in  the  Northern  Hemisphere.   So  in  almost  all  cases 
the  3F/3AC  values  for  the  monthly  data  are  lower  than  for  the  annual  data. 

The  above  discussion  shows  that  under  certain  conditions  dF/dTs  is  approxi- 
mately equal  to  3F/3TS  whereas  for  other  conditions  they  differ  appreciably.   Since 
in  this  study,  the  only  variables  affecting  F  are  assumed  to  be  T  and  A  ,  it  can 
be  shown  that  dF/dTs  should  differ  from  3F/3T   except  under  two  conditions  when 
either  3F/3AC  or  dAc/dTs  is  equal  to  zero. 
Mathematically  this  follows  from 

3F         3F  . 
dF  =  (— )  dT„  +  (— )  dk  (3) 

s  c 

or  dF  =  _3_F  +  d¥_     dA,  /4) 

dTs   3TS  (  3A^dT^ 

From  (4)  dF  =  3F    ,    .     ... 

-j=       -r=     only  when  either 
dTs   3TS 

_3_F  =  0,  or  dA  =  0 
c  s 
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But  our  calculation  has  shown  that  9F/^A  is  a  significantly  large  number  (in 
absolute  magnitude);  therefore  for  dF/dT  to  be  equal  or  very  close  to  S>F/dT  , 

o  S 

dA  /dT  must  be  equal  or  very  close  to  zero,  i.e,  A  and  T  must  not  be 
correlated.   This,  indeed,  is  the  case  except  in  the  annual  zonal  case  for 
both  Northern  Hemisphere  and  ocean  where  high  correlation  was  found  to  exist 
between  A  and  T   (3^=0.84  and  0.76  respectively).   And  it  is  exactly  these  two 
instances  in  which  dF/dT  and  9F/9T  are  significantly  different.   We, 

S  o 

however,  hesitate  to  attach  any  physical  significance  to  this  finding  since 
the  sample  size,  only  9,  is  small. 

In  order  to  examine  in  more  detail  the  sensitivity  that  clouds 
demonstrate,  Table  4  was  generated.   All  data  employed  to  produce  numbers  in 
Table  4  are  identical  to  those  used  for  Table  1  except  that  the  cloud  fraction 
is  based  on  London  (1957).   London's  cloud  amounts  are  lower  at  all  latitudes 
compared  to  those  of  Berlyand  et  al.  (1980)  used  for  Table  1.   The  differences 
in  c  between  Table  4  and  Table  1  are  very  significant  especially  for  the 
Northern  Hemisphere  where  differences  in  b=(9F/9Tg)  are  also  very 
noticeable.   This  demonstrates  again  the  delicate  nature  of  cloud  sensitivity 
to  the  cloud  data  set  used,  as  was  discussed  earlier. 

An  alternative  way  to  estimate  b  and  c  is  to  look  at  the  variation  of 
longwave  radiation  with  respect  to  surface  temperature  and  cloud  cover  at  each 
latitude  zone.   Table  5  lists  these  values,  using  surface  temperatures  from 
Warren  and  Schneider  (1979),  NOAA's  SR  longwave  fluxes  from  Gruber  and  Winston 

(1978)  and  fractional  cloud  data  from  Berlyand  et  al.  (1980).   To  ensure  the 

2 
reliability  of  these  estimates,  only  those  values  with  r  >0.85  are  kept. 

Because  of  the  small  annual  variation  of  temperature  at  low  latitudes  and 
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uncertainties  in  the  data,  which  lead  to  low  correlations,  it  is  primarily 
these  latitudes  that  are  excluded  by  the  above  constraint.   Further,  if 
the   standard   error  of  a   coefficient   (SE)   is  comparable  to   the 
magnitude  of  the  coefficient  itself,  the  true  value  of  such  a  coefficient 
must  be  viewed  as  unreliable.   It  is  quite  clear  that  b=(dF/dTs)  is  not 
constant  but  has  significant  meridional  variations,  generally  decreasing  from 
equator  to  pole.   Similar  variations  of  b  with  latitude  occur  in  the  Southern 
Hemisphere.   Such  variations  have  also  been  found  by  Warren  and  Schneider 
(1979),  who  used  the  Ellis  and  Vonder  Haar  (1976)  radiation  data  set,  and  by 
Ramanathan  (1977)  from  simulations.   However,  the  larger  value  of  b  obtained 
at  higher  latitudes  in  the  Southern  Hemisphere  by  Warren  and  Schneider  (1979) 
is  contrary  to  our  result. 

A  decrease  of  b  toward  the  poles  might  be  expected  from  the  meridional 
variation  of  the  effective  emission  temperature  (T  )  of  the  Earth-atmosphere 
system.   In  general,  Te  decreases  with  latitude,  at  least  for  latitudes 
greater  than  20-30  .   If  the  outgoing  longwave  flux  is  represented  by  (fT   , 
where  U  is  the  Stefan-Bol tzmann  constant,  then  dF/dT  would  decrease  with 
increasing  latitude.   Since  the  effective  emission  temperature  is  positively 
correlated  with  surface  temperature,  one  might  also  expect  that  b  =  dF/dT 
would  also  decrease  toward  the  poles. 

3-2   Relationship  between  surface  and  planetary  albedo  for  clear  skies. 

Regression  equations  linking  the  planetary  to  surface  albedo  for  clear 
skies  were  derived  by  employing  data  obtained  from  model  calculations  as 
described  in  Section  2.   These  relationships  were  then  tested  with  real  data 
from  satellite  observations  and  ground  measurements.   The  method  of  assessing 
the  derived  equations  was  to  compute  the  root-mean-square-error   (RMSE)  of  the 
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planetary  albedo  defined  by, 


RMSE  = 

A4 

where  n  is  the  number  of  observations ,  oC  >  t-he  observed   planetary  albedo, 
and  o^the  predicted   planetary  albedo  defined  as  follows, 

c£   -  0.0587  +  0.7301  JT  (6) 

tf(  =  0.0410  +  0.7610  If  (7) 

where  f   denotes  the  surface  albedo.   It  is  necessary  to  point  out  here  that 
the  (G\.   \    )  relationship  in  Eq .  (6)  was  derived  for  all  solar  zenith  angles, 
whereas  Eq .  (7)  was  only  for  a  solar  zenith  angle  of  40°. 

Both  Eqs.  (6)  and  (7)  show  that  at  low  surface  albedos  the  planetary 
albedo  is  greater  than  the  surface  albedo;  this  is  a  result  of  the  importance 
of  the  contribution  of  Rayleigh  scattering  to  the  planetary  albedo  when  the 
surface  albedo  is  small.   At  large  surface  albedos,  the  planetary  albedo  is 
less  than  the  surface  albedo;  this  is  due  to  the  depletion  of  the  solar 
radiation  by  atmospheric  absorption,  which  reduces  the  reflected  solar 
radiation  reaching  the  top  of  the  atmosphere.   There  is  a  crossover  point  at 
which  the  planetary  and  surface  albedos  are  equal.   For  Eq.  (6),  this  occurs 
at  an  albedo  value  of  0.22,  for  Eq.(7),  at  an  albedo  of  0.17. 

Table  6  shows  that  the  RMSE  is  0.05,  except  when  the  Hummel  and  Reck 
(1979)  values  of  |f  are  used  in  Eq.(7),  in  which  case  RMSE  is  0.06.   This  may 
be  attributable  in  part  to  the  fact  that  Hummel  and  Reek's  seasonal  data  were 
interpreted  as  representing  the  mid-seasonal  month. 

Table  7  lists  two  regression  equations  derived  from  the  observational 
data,  using  respectively,  Robock  and  Hummel  and  Reck  for  surface  albedo.   In 
each  case  the  standard  error  is  0.05.   These  two  equations  are  to  be  compared 
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with  Eqs.  (6)  and  (7)  which  were  previously  derived  using  theoretical  data. 
Notice  that  both  the  intercept  and  slope  of  the  first  equation  in  Table  7  are 
comparable  to  their  counterparts  in  Eq  (6).   This  indicates  again  that 
Robock's  data  are  in  somewhat  better  agreement  with  the  theoretical  equations. 

Fig.  1  shows  the  plot  of  the  four  individual  regression  lines,  two  from 
the  simulations  and  two  from  the  observations.   At  the  higher  and  lower  ends 
of  both  oCand  f     ,  there  is  a  difference  of  about  0.05  between  the  highest  and 
lowest  of  these  lines.   In  the  middle  range,  (.0=   0.2  to  0.4),  the  differences 
among  the  regression  lines  are  very  small.   Overall,  the  line  representing  Eq. 
(6)  agrees  best  with  the  other  three  lines.   Further,  these  three  lines 
certainly  would  fall  inside  the  domain  bounded  by  upper  and  lower  lines  drawn 
parallel  to  the  line  representing  Eq.  (6)  to  distance  equal  to  ±  0.05. 

A  somewhat  similar  study  was  made  by  Preuss  and  Geleyn  (1980).   For  given 
model  atmospheres,  and  using  the  two  stream  approximation,  they  calculated  the 
reflected  solar  radiation  for  two  different  surface  albedo  values.   For  each 
model  atmosphere  they  then  derived  a  linear  relation  between  the  planetary 
albedo  and  surface  albedo.   The  average  values  of  the  intercept  and  slope  of 
their  equation  are,  respectively,  0.101  and  0.59  for  June  and  0.108  and  0.57 
for  January/February.   These  numbers  differ  quite  significantly  from  our 
observational  and  theoretical  results  discussed  above. 

The  agreement  of  the  NOAA  SR  albedo  observations  with  the  theoretical 
regressions  suggests  that  although  the  SR  albedos  are  based  upon  measurements 
in  the  visible  region  of  the  spectrum  they  yield,  at  least  for  zonal  averages, 
reasonable  values  for  the  total  albedo. 

Recently,  in  connection  with  climate  modeling  studies  some  interest  has 
been  generated  in  the  magnitude  of  the  surface  albedo  of  the  Sahara  desert 

-9- 


region.   Based  on  observations  of  16  target  areas  (^^500k.m  x  500km  each)  over 
the  Sahara  desert  area  for  a  period  of  12  days,  the  ERB  (NIMBUS  7)  gave  an 
average  minimum  albedo  of  36.7%.   If  one  inserts  this  value  into  Eq.  (6),  one 
obtains  a  surface  albedo  of  around  0.42.   Similarly,  if  one  uses  the 
regression  of  ^fonoC,  i.e.,y=  -0.0724  +  1.34640L,  one  also  obtains  a  surface 
albedo  of  0.42,  since  in  this  case  qJC and  y fell  nearly  precisely  on  the  same 
straight  line.   Since  the  NIMBUS  7  observations  are  close  to  local  noon  and 
since  there  is  a  tendency  for  the  surface  albedo  to  increase  with  increasing 
solar  zenith  angle,  the  actual  mean  daily  surface  albedo  is  probably  a  few 
hundreths  higher  than  0.42.   Interestingly,  an  even  higher  value  is  predicted 
by  Preuss  and  Geleyn  (1980)  equation.   These  albedos  are  much  higher  than 
those  usually  attributed  to  desert  regions  (e.g.,  Budyko,  1974,  p. 55,  gives  a 
value  of  0.28,  and  Hummel  and  Reck,  1979,  give  values  of  0.25  to  0.31  for  the 
Sahara  region  ),  although  Otterman  and  Fraser  (1976)  have  derived  a  value  of 
0.44  from  Landsat  observations. 

4.   Conclusions 

The  results  presented  here  show  that  the  cloud  sensitivity  coefficient  ( 
5F/3A  )  as  determined  from  multiple  regression  analysis  of  highly  averaged 
data,  is  strongly  dependent  on  the  radiation  and  cloud  data  sets  used  in  the 
analysis  and  also  on  whether  one  uses  annual  or  monthly  means.   To  a  much 
smaller  extent  this  is  also  true  of  the  temperature  sensitivity  coefficient. 
There  are  a  number  of  implications  for  climate  models  and  climate  theory. 

1)   Presently  available  radiation  and  cloud  data  sets  yield  different 
values  for  longwave  radiation  sensitivities  causing  uncertainties  in  the  true 
values  of  these  coefficients. 
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2)  Because  of  the  uncertainty  in  the  value  of  dF/dT  to  be  used  in 
energy  balance  climate  models,  the  results  of  such  models  are  subject  to 
uncertainty  (see  Warren  and  Schneider,  1979  for  sensitivity  of  climate  model 
results  to  dF/dT  ).   Our  results  suggest  that  rather  than  using  a  single  value 
of  the  coefficient  dF/dT  in  such  models,  a  latitude  dependent  or,  even 
better,  a  temperature  dependent  coefficient  should  be  used. 

3)  The  strong  sensitivity  of  9F/9A  to  radiation  and  cloud  data  sets  and 
to  time  averaging  period  suggest  that  it  was  fortuitous  that  Cess  (1976) 
obtained  a  value  of^-"  -90  W"m    for  ^F/3A  —  a  value  almost  exactly  balancing 
the  cloud  albedo  effect.   It  would  appear  that,  if  one  were  to  try  to  obtain  a 
mean  hemispheric  or  global  value  of  9F/9Ac  from  a  regression  of  the  sort  used 
here,  it  would  be  more  appropriate  to  use  the  monthly  data,  which  allow  for 
seasonal  and  latitudinal  variations,  rather  than  the  annual  data,  which  allow 

only  for  latitudinal  variations.   In  almost  all  cases  the  3F/9A  values  for 

J  c 

the  monthly  data  are  lower  than  for  the  annual  data  and  are  well  below  the 
value  required  to  balance  the  cloud  albedo  effect. 

4)  To  obtain  better  estimates  of  these  important  sensitivity  parameters, 
more  accurate  earth  radiation  budget  and  cloud  climatological  data  are 
required. 

A  simple  linear  relationship  has  been  derived  between  the  planetary 
albedo  and  the  surface  albedo  for  the  case  of  clear  skies.   Such  a 
relationship  has  important  applications  in  climate  modeling  and  in  the  use  of 
satellite  observations  for  cloudiness  estimates  and  surface  energy  budget 
determinations.   Application  of  this  relationship  to  satellite  albedo 
observations  over  the  Sahara  desert  region  yields  a  surface  albedo  of  0.42, 
significantly  higher  than  previous  estimates. 
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Table  6.   Observational  data  sets  used  to  test  the  derived  regression 
equations^  O^and  IT are  respectively  the  planetary  and  surface  albedo. 


Equation 


RMSE 


Data  Sets 


<?C  =  0.0587  +  0.7301  ^ 
<K  =  0.0410  +  0.7610  V 


0.05 
0.05 


NOAA  SR  («(.),  Robock  (If) 
NOAA  SR  (4C),  Robock  (flf) 


(7<v  =  0.0587  +  0.7301^ 
^  =  0.0410  +  0.7610  J' 


0.05 
0.06 


NOAA  SR  CsO,  Hummel  &  Reck  (JO 
NOAA  SR  (oO,  Hummel  &  Reck  QT) 


Table  7.   Regression  equations  derived  from  observations. 


Equation 


Standard  Error 


Data    Sets 


<£  =  0.0597  +  0.7560  ^ 
*i=   0.0784  +  0.6760  % 


0.05 
0.05 


NOAA  SR  (©(.),  Robock  (jf) 

NOAA  SR  (*C),  Hummel  &  Reck  if) 
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Figure  1 


Plot  of  four  regression  lines,  two  each  from  simulations  and 
observations.   The  thin  dashed  and  solid  lines  show  respectively 
the  simulations  for  all  solar  zenith  angles  and  for  a  solar 
zenith  angle  of  kO      only.   The  heavy  dashed  line  is  from  observa- 
tions of  NOAA-SR  and  Robock  data.   The  dashed-dotted  line 
delineates  observations  using  NOAA-SR  and,  Hummel  and  Reck  data. 

AB  represents  the  line  of  o(.  =   q    >    when  there  is  no  atmosphere. 
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